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ABSTRACT 

An    investigation    of   conventional    and   broken -symmetry    scaling, 
and   their    effect    on   the  ratio   of   the  neutron-proton    electric    form 

factors,    G^   /G„    ,    was   carried  out   using   measurements    from   elastic 

En      Ep 

electr on-deuter on   and   electron-proton    scattering    in   the  range   of 

2  2  -2 

momentum  transfers,    q    ,    of  O.lO^q  ^0.80   F      ,    and    electron    scat- 
tering  angles    of   45      to   120    .      Within    experimental    errors    of 
available  data    in    this    range   of  momentum   transfers    and  angles,    it 
was    concluded   that   both  forms    of   scaling   as   applied   to  the   deuteron 

form   factor    have   negigible    effect    on    the  ratio  G„   /G^ ,    and   the 

En      Ep 

application    of    either    scaling    law  to   the  proton    form   factor    results 
in   nearly    identical   corrections    to   this    ratio,    thus   allowing   no 

distinction  between   them,    and   the  agreement    of   the   experimental 

2 
data   with   the   neutron-electron    interaction   slope  at    q     —  0    is    not 

influenced  by   the   choice   of   scaling    lav/.      A  conclusive   test    for 

the   scaling    laws    would  require  an    improvement    in    experimental 

accuracy    by   at    least   an    order    of  magnitude,    a    technically    nearly 

impossible  requirement. 


TABLE  OF   CONTENTS 

I.  INTRODUCTION 6 

II.  THEORY 10 

A.  GENERAL 10 

B.  SCALING  APPLIED  TO  MAGNETIC   CORRECTIONS l4 

C.  DERIVATION   OF    RELATION    FOR    G^   /G^_    WITHOUT 

En      Ep 

SCALING- 17 

III.  REEVALUATION   OF   EXPERIMENTAL    DATA 23 

A.  DATA 23 

B.  ROSENBLUTH  PLOTS 23 

C.  MAGNETIC  CORRECTIONS 24 

D.  CALCULATION  OF  G^  /G^   WITHOUT  SCALING 24 

En   Ep 

E.  CALCULATION  OF  G,,  /G^_   WITH  SCHUMACHER  SCALING--25 

En   Ep 

IV.  PRESENTATION  OF  RESULTS,  DISCUSSION,  AND  CONCLUSIONS-26 

A.  PRESENTATION  OF  RESULTS  AND  DISCUSSION 26 

B.  CONCLUSIONS  28 

APPENDIX  A.  TABLES  30 

APPENDIX  B.  FIGURES  37 

APPENDIX  C.  STATISTICAL  TREATMENT  OF  DATA 43 

BIBLIOGRAPHY -- 46 

INITIAL  DISTRIBUTION  LIST 47 

FORM  DD  1473 48 


LIST   OF   TABLES 

I.  Exper  imental    Data 30 

2 

II.  Mott  Cross  Sections,  G  ,  Theoretical  Cross  Sections 

P 

and  Normalization   Factors 31 

M  E 

III.  D    ,    D     and    D     Computed  from   Feshbach-Lomon   Wave 

c    M       M 

Functions 32 

IV.  Deuteron  Magnetic  and  Electric  Form  Factors  from 
Rosenbluth  Plots 33 

V.  Effect  of  G    on  G    Based  on  Rosenbluth  Plots 34 

VI.  Magnetic  Corrections  and  Effect  on  G   /G    as 

En      Ep 

Computed  by   Stewart 35 

VII.  G^   /G_      Computed  Without   Scaling,    with  Proton 

En   Ep 

Scaling,  with  Schumacher  Scaling,  and  as  Computed 

by  Stewart 36 


LIST  OF  FIGURES 

1.  Rosenbluth  Plot  for  Deuteron 37 

2.  Electric  Form  Factors  for  Deuteron 38 

3.  Magnetic  Form  Factors  for  Deuteron 39 

4.  Ratio  of  Magnetic  to  Electric  Form  Factors  for 
Deuteron 40 

5.  Comparison  of  G^   /G^_  Computed  Without  Scaling, 

En   Ep 

and  with  Proton  and  Schumacher  Scaling ■ 4l 

6.  Comparison  of  G   /G„   Computed  Without  Scaling  and 

G^  /G^_  as  Computed  by  Stewart 42 

En   Ep        ^       ' 


I.   INTRODUCTION 

The  elastic  scattering  of  electrons  by  protons  is  described  in 
the  first  Born  approximation  by  the  Rosenbluth  equation  [l].   In 
this  approximation  a  single  virtual  photon  is  exchanged  between  the 
electron  and  the  proton  and  the  structure  of  the  proton  is  repre- 
sented by  two  electromagnetic  form  factors  which  are  functions  of 

2 

q    ,    the   square   of   the   four-momentum   transferred  by    the   virtual 

photon   [2].      Janssens,    et.al.,    L3]    made   absolute  measurements    of 

the   elastic    electron-proton    cross    sections    with   an   accuracy   of 

2  -2  -15 

about    4%   in   the    q      range  from   4.0   to  30.0   F      (1   F   =    1   Fermi    =    10        m) 

These   measurements    verified   the   Rosenbluth   equation    for    values    of 

2  -2 

q     up    to   22.0   F      ,    and   determined  the   form   factors    of    the  proton 

to  a   precision    exceeding    that    obtained   in  previous    experiments.      It 

was    also   found   that    the    charge  and  magnetic    form   factors    of   the 

2 
proton   have   the   same    dependence   on    q     within    the  accuracy   with 

which   they   were   determined. 

Hundreds    of   other    elastic    electron-nucleon    scattering  measure- 

2 
merits    are  also  available,    extending    over    a   very  wide  range   of    q 

and   angle   L4] .      Empirically,    the   data    seem   to   exhibit   remarkable 

simplicity,    the   electric   and   magnetic    form  factors    of   nucleons 

obeying    the    "scaling    law"; 

2  GMo(q2)         GMn(q2)  2 

%<<>>* 'IT =   °SL;    GEn<q    >    =    °         ™ 

P  n 

2 
with  the  q   dependence  being  well  described  by  the  "dipole  fit"; 


GSL  ~  GDipole   -  T~2  (2) 

(1  +  qVM^) 

2  -2 

where  M   =  90.25  F   •   It  has  long  been  clear  that  the  electron- 
D 

proton  data  show  systematic  oscillations  about  the  dipole  fit 
[5}°]»  but  the  amplitude  seems  to  vary  with  laboratory  and  (2)  is 
regarded  as  a  remarkably  good  phenomenological  fit. 

Following  the  same  general  procedure  as  Rosenbluth  did  with 
the  proton,  Jankus  [7]  derived  the  cross  section  for  the  elastic 
electron-deut er on  (e-d)  scattering  from  a  point  deuteron  in  the 
first  Born  approximation.   His  results  differ  qualitatively  from 
the  proton  cross  section  in  one  important  way,  i.e.,  there  are 
three  form  factors  characterizing  the  deuteron  while  two  suffice 
for  the  proton.   This  is  in  agreement  with  the*  work  of  Rlaspr  and 
Jaksic  [8J.   Their  study  showed  that  the  cross  section  for  the 
scattering  of  a  relativistic  electron  from  a  potential  with  spin  J 
contains  2J  +  1  form  factors.   The  three  deuteron  form  factors 
correspond  to  the  charge,  quadrupole,  and  magnetic  moments  in  the 

2  r       1 

static    limit    q     ~*   0 .       Dr  icker    and   Hand   L  9  J    experimentally    verified 
for    the   deuteron   an    emperical    relation    between    the   form   factors 
which    is    analogous    to    that    for    the  proton.      The    scaling    law   for    the 
deuteron   may    be   written; 

G  M 

GED~5T    .    !>{,-    1-71    =5^^      .  (3) 

D  p 

Based  on  these  scaling  laws,  Stewart  [10]  in  a  study  of  the 

neutron  charge  from  factors  at  very  low  momentum  transfers,  found 


that: 

(a)  The   Peshbach-Lomon   wave   functions    [ll],    together    with 

relativistic   corrections,    removed  an   apparant    discrepancy  between 

2 
the  neutron-electron    interaction    slope  at    q     —   0,    and    the   slope 

given   from    values    of   G         (obtained  by   electron    scattering)    in 

the  range   0.10  ^    q      ^    0 . 80   F~     . 

(b)  Within  the  relatively  large  errors  propagated  into  G   ,  the 
Partovi  wave  functions  L12j  with  relativistic  corrections  applied, 
are  in  disagreement  with  the  neutron- electron  interaction  slope. 

The  important  conclusion  from  Stewart's  experimental  results 
is  that  the  neutron  has  a  non-zero  charge  form  factor,  implying 
a  charge  distribution  within  the  neutron. 

Schumacher  L4J  however,  points  out  that  the  theoretical  sit- 
uation is  not  satisfactory.   The  scaling  law  (1)  led  to  the  in- 
vention [13]  of  the  phenomenlogical  symmetry  now  called  SU(6)  , 
but  a  deeper  interpretation  appears  to  be  lacking.   Likewise,  the 
dipole  formula  (2)  corresponds  to  a  double  pole  of  mass  843  MeV 
instead  of  a  combination  of  single  poles  at  masses  of  realistic 
vector  mesons.   Although  it  was  hoped  that  the  oscillations  around 
(2)  might  be  explainable  in  dynamical  terms,  Goitein,  Dunning  and 
Wilson  [5]  showed  that  existing  theories  were  not  adequate  to  do 

so  and  that  form  factor  dominance  by  vector  mesons  led  to  trouble 

2  .  . 

at  high  q  .   There  has  been  little  change  in  the  situation  they 

described,  e.g.,  the  best  subsequent  models  still  not  being  for- 
mulated such  that  an  isospin  reflection  connects  realistic  de- 
scriptions of  the  proton  and  neutron  [4].   Because  of  this  awk- 
ward situation  regarding  suitable  dynamical  models  for  the  form 


factors,  it  is  desirable  to  ask  the  data  to  provide  as  many 
answers  as  possible  to  questions  which  can  be  formulated  without 
a  committment  to  the  restrictive  details  of  particular  models. 

Schumacher  examines  the  symmetry  content  of  (1)  and  its 
dynamical  implications,  and  finds  that  the  form  of  the  deviations 
from  (1)  suggests  a  particular  type  of  symmetry  breaking,  well 
known  in  the  physics  of  strong  interactions.   From  these  consider- 
ations a  broken  symmetry  theory  of  the  form  factors  is  formulated, 
Using  the  world's  e-p  data,  Schumacher  obtains  several  measures 

of  symmetry  breaking,  which  imply  an  oscillating  G    whose  sign 

En 

and  magnitude  agree   with    the  Krohn-Ringo   slope    Ll43.      He    further 
derives    a    set    of   scaling    laws    which   relates    the   form   factors    G 


Mp 


and  G..      to  G_     and  G_    .      These   relations    are 
Mn  Ep  En 


G         —  LL    G         +  u-    G 
~Mp         ^p^Ep  n"En 

(4) 

Gw      =110,,      +  p,    G„ 
Mn  n    Ep  p    En 

where  G^    =   G^_    +   G„„    is    the    "symmetry-conserving"   part   and 
Ep  ES  EV  r 

G^      =   G„_    =   G„„    is    the    "symmetry-violating"   part    of   the   form 
En  ES  EV  r 

factors.      These   conclusions    are  primarily   model- independent    and 
do   not   require   pole   dominance   of   the    spectral    functions. 
It    is    then,    the  purpose    of   this    thesis    to: 

(a)  Test    the   two   forms    of    the   scaling    laws    using    the    data    of 
Bumiller    and   others    [l9*]    in    the   range   0.10  <■    q      <■    0.80   F~    ,    and 

(b)  To    introduce   a   method   of   finding    the   ratio  G_    /G_      without 

En   Ep 

using  the  scaling  laws,  thus  finding  information  about  G„  . 

En 


II.   THEORY 

A .   GENERAL 

Jankus  [7],  found  that  for  the  deuteron 


dcr 
da 


V^'Mott 


L  2E       .    2    0 

!   +  —  sm     - 

>  MDC 


-ir 


+   I  T?(l    +    2(l+r?)tan2   %)    D^(q2) 


2'       W 


2.     2X       8      2    2,    2, 
Dc(q    )+    9   1?     DQ(q    ) 


(5) 


wher  e 


Dc(q2)  =I   [u2(r)+  w2<r) 


J0(T)dr 


(6) 


is    the   charge   structure   factor 


,    2X         6  NTT     2 

Vq    )    =         2      % 

q  ~o 


u(r)w(r ) 


w   (r) 


J2(T)dr  (7) 


is    the   quadrupole   moment    structure    factor ,    and 


°>2>=I 


u    ( r  )  -%  w   ( r  ) 


j    (T)dr    +  —  \ 


u(r )w(r) 


2/     ^ 

,w    (r) 


52(T)at 


(8) 


accounts  for  the  contribution  of  the  intrinsic  magnetic  moments 
of  the  proton  and  the  neutron  to  the  scattering  process,  and 


D^(q2H(   w2(r) 


J0(T)  +  J2(T) 


dr 


(9) 


is  the  magnetic  contribution  to  the  scattering  process  arising 
from  the  convection  of  charge  in  the  deutron.   Together  (8)  and 


10 


(9)    are 


DM(q    )    =% 


(U    +   [i     )2DM    +    D 
v    p         n'       M  M 


The  static    limits    (q  -•   0)    are 


Dc(q2) 


vq2) 


Vq    > 


2^ 

M    Q 
D 


1    "iPD 


(10) 


2       D 


The   j      and   j      are   spherical    Bessel    functions,    T    =    |q    |r/2,    Q    is 
the    quadrupole   moment,    P      the  percentage   D  state,    and   the    u(r)    and 
w(r)    are   suitably    chosen   wave   functions    of   the    deuteron. 

Adler    Ll5Jj    shows    that    in    the   impulse   approximation,    the 
charge,    quadrupole,    and   magnetic   form   factors    are   given   by 
(along   with   their    static    limits): 

GC   3    <GEp    +    GEn>     DC  2— T1  (11) 


and 


V"    (%    +   GEn>     DQ 


"MDQ 


(12) 


and 


M 


E  M 

G=(G         +   G       )    D      +    (G         +G      )2D    - 

M         v     Ep  En'       M  v    Mp  Mn ;       M 


n 


"UD      U3) 


Ep  En '       M         v  >lp 

F 

where  the  D's  are  given  in  equations  (6)  through  (9).  The  step 
from  a  point  deuteron  to  a  finite  size  deuteron  is  thus  carried 
out  by  the  introduction  of  the  free  nucleon  form  factors.  This 
step    is    justified   by    the   impulse   approximation. 
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Rewrite   equation    (5)    as 


where 


/d7\    ,/dCr  ^  2         A/    2.         nt    2,  2   9  ,  _    . 

v  v        /  Mott 


,    2,  2         8      2    2         2  2 

A(q  )    =Gc  +  -r7  GQ  +  -r?  GM 


and 


T? 

2 

-  q 

,    2    2 

4mdc 

B(q2) 

■    ^  >r>    r2 

~   3^  GM 

(5c) 


where   the   approximation    1    +  rj    —   1    has    been    used.      Equation    (13) 
can   be   simplified  considerably   by    scaling,    i.e., 


and 


G         =   U  'G       ,      u.  '    =    1  .71 
MD  D    CD'  D  ' 


G  G 

Mp    _      Mn    _ 

U-p  l^n  Ep' 


For  the  case  of  the  deuteron 


G        =    (G      +   G      ){i,  '    D    .  ( 14) 

MD         v    Ep         En^D      C  v        ' 


Then 


or 


2  ,  s2r    2/n       2  ,2      4  ,2         2    9N       8      2    2    t 

GD  =    (GEP+  GeJ    Cdc<1  +  3^  +   3^0  tan     2>  +  9.^    °Q  ] 

4=    <GEP+    GEn)2    ^  <l5> 

2 

where   F'       is    called   the   deuteron    structure   factor    and    is    the   same 

as    the  result    of   Jankus    for   point    nucleons    in    the    deuteron.      From 
(15) 

.  ^=p.(x+^\.  (16, 

Ep  \  Ep  / 

Again    consider    equation    (5)-       If    the    quadrupole   term    is 
neglected    (a    contribution    of   order    10         or    less)    and    1    +  T\   ''      1 


12 


is    used,    (5)    becomes 


/dtf\     .  /dCT  \  2         2         2  _  2    0X    2 

.2         2         2.        4      ^2    A       2    0  ,,„. 

=    (GE    +  ^V*   FGM    tan      2     •  (17) 

^ d  a  a  /  /d  a  "\          *».  ^-  * 

quation,    using     /  — — j   //  — )               as    the  ordinate 

^  Mott 


From   this    e 


2 
and   tan      9/2   as    the  abcissa,    a    linear    plot    can   be   made    ( "Rosenbluth 

2  2 

plot").      Then   G        is    given   by    the   slope,    and  G    ,    can    be   found  from 

the    intercept.      This    method,    however,    presumes    that   absolute  mea- 
surements   of    the   deuteron    cross    section  are   available.       Since 
Stewarts    measurements    are   relative,    available    data   must    be  an- 
alyzed  in   an    attempt    to   make   the    experimental    cross    sections 
absolute    in    a    first    order    sense. 

This    task    can    be  accomplished  by    considering    the  work   of 
de   Vries    [l6]    and    Buchanan    [4] .      De   Vr ies    collected  all   available 

proton    data,    and   made  a    fit    to   these    data    using   adjustable   para- 

2 
meters    so   that   a    minimum  was    obtained    in   a  X      fit.       Buchanan 

collected   the   world's   proton,    neutron,    and   deuteron    data,    and 

updated   them.      He    checked   the   de   Vries    b'    fit    using    the   reanalyzed 

data,    and  concluded   that    this    fit    is    the   best    to    describe    the 

proton    form   factors. 

Stewart    [lo]    measured 

-0fr)7'(a-)r- 

2 
If  it  is  assumed  that  the  de  Vries  b'  fit  gives  G   absolutely, 

P 

then 

(-)P  =(vabs(-p    •  <15'  , 
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If   the   ratio 

.  ■  (C  («■)  7 

is    formed,    N   can    be   considered    to  be   the    factor    required   to   nor- 
malize  the   experimental    data    to  absolute    data    in   a    first    order 
sense.      Then,    since    the   ratio   R    was    measured 


do\theor       /d0_]exp 


(N)  (21) 


can  be  considered  to  be  absolute  in  the  first  order.   The  sub- 
script theoretical  is  used  to  distinguish  absolute  in  the  first 

order  from  absolute. 

2 
Thus  Rosenbluth  plots  for  the  deuteron  may  be  made,  and  G 

ED 

2  .  .   . 

and  G    can  be  extracted.   If  the  scaling  law  holds,  then  within 

MD  * 

experimental    error 

GMD 


M-  'G 
D  ED 


B.   SCALING  APPLIED  TO  MAGNETIC  CORRECTIONS 

2       2 
To  eliminate  the  magnetic  contribution  to  G   and  G  ,  the 

magnetic  contribution  is  given  by  equation  (15)  as 

CD  =  f^D2'1  +  2  ta"2!>  <22> 

so  that  2 

2     Gp  exP 

GED  "  (1+CD)  (23) 

with 

where  the  quadrupole  terms  has  been  neglected. 
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For  the  proton,  the  magnetic  contribution  is  found  by  con- 


sidering                                        2  2 

G2    =   rg-   +  T7  -53EL-    (i    +   2    tan2   §)  .             (24) 

p         1+77               p  l+n       V                                2'                v    ^> 

P  P 


Scaling  assumes  G.^  =  M-   G,,  ,  where  u-   =  2.973,  so, 

Mp  p      Ep  p 

gp  =  i^- <x  +Vp(1  +  2  tan2!)}-  (24a) 

The   term 

C      =  T)   |i       (1    +   2    tanZ  |)  (25) 

p  'p   p    x  2'  v       ' 

is  then  the  magnetic  correction  term  for  the  proton.   Thus 

2  Gd(1    +  V 

G2     =  -2= E-  (26) 

Ep  1    +   C  v       ; 

P 

where  z        \  /        \ .,    .. 

<--(*r'(«r- 

From   his    experiment,    Stewart    LIOJ    calculated 

GEn    _    1     GED 

G  F      G 

Ep  D      Ep 


and 

eva 

luat 

es    f r  om 

his    data 

<4 

E 

KD 

4, 

F 

K 
P 

(27) 


where   the   E  and  F   contain   radiative    corrections,    Mott   cross 
sections,    target   parameters   and  experimental    counting   rates. 
Of    interest    here    is    the  ratio  K^/K    ,    which   represents    the 
magnetic    correction    terms,    it    is    given   by 

Kn  1   +   C 

K  (1+C   )(1+T)    )  K       ' 

p  x  D  p 


15 


where   C     and   C      were   defined   by    (22)    and    (25).      This    ratio   contains 
the   scaling    law  assumptions. 

Next    consider    the    scaling    law  proposed  by    Schumacher    C4J . 
Specifically 


G         =   U      G         +   u      G 
Mp  p      Ep  n      En 


(29) 


G         =  \i       G         +  II      G 
Mn  n      Ep  p      En    . 

For    the    deuteron,    consider    similiar    scaling    laws 


G        =  u  f    G        +  u      G 
MD        PD      ED        Pn      En 


(30) 


G         =   U      G         +   n  '    G 
Mn  n      ED  D      En    . 


Then  equation  (17)  becomes,  neglecting  terms  containing  G   , 


(G^    +   G^)2    D2Q    [l 


.2 
'ED 


Ep  En 


2  ,2 

3^D 


2    0 
2 


+  ^A^(1  +  2tan2|K| 

Ep      En    V"  C  /J 


(31) 
(32) 


Now   define 


cd  =  (!vd  + 


4   GEP   GEn     »i»u,W 
3   GRD+GRn         DC        D/  V 


1+2    tan 


2    0 


(32) 


Ep      En  C 

Similar  ily,    for    the   proton    equation    {2L\)    becomes,    with    the   use 


of   equation    (29), 
.2 


.2         G_Ep_ 


G      = 


P       1-^7, 


I   H2+   2t)   a   ii     -^    )/l    +   2    tan2   | 

p  p      p  p  nGEp  A  2 


(33) 


and   define 


(Vp 


P        I  P     P  P   P    n    G£p 


2  e\ 

n      -j      . 


(34) 


Thus    when    using    Schumacher's    scaling,    the   magnetic    correction 


ratio  becomes 


<1  »  c;> 


(1+Cp(l4tjp      ) 


(35) 
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It  should  be  pointed  out  that  equations  (32),  (34),  and  thus  (35) 

contain  the  form  factors  G„   and  G^  ,  both  of  which  are  sought 

En       Ep  v 

from  experiment.   To  use  the  readily  available  values  from  the 

de  Vries  b*  fit  would  be  combining  the  dipole  fit  with  the  broken 

symmetry  theories,  however  these  are  needed  for  comparison  of  the 

scaling  laws.   If  the  de  Vries  b'  fit  values  of  G„   and  G„   are 

Ep       En 

2 
incorrect  at  low  q  ,  their  influence  on  equation  (35)  is  quite 

small,  so  they  will  be  used  here. 

A  comparison  of  equations  (28)  and  (35)  gives  a  measure  of  the 
effect  of  the  two  scaling  laws  on  the  data  of  Stewart.   These  re- 
sults are  shown  in  Table  VI. 


C.   DERIVATION  OF  AN  EQUATION  FOR  G^  /C^      WITHOUT  SCALING 

En   Ep 

This  section  has  a  threefold  objective: 

(a)  Derivation  of  an  expression  for  G  /G   without  the  use 

v     '  '  En      Ep 

of    scaling    laws, 

(b)  Testing   and   comparison    of   the    values    of   G_    /G„      obtained 
x     '  =>  r-  En      Ep 

2  -2 

from    the    equation    derived    in    (a)    above,    for    0.1   ^    q      ^    0.8    F      , 

using    de   Vries    b'    fit,    or    proton   scaling   wherever    required. 

(c)  Testing   the    effect    of    Schumacher's    scaling    on   G^    /G__ . 

En   Ep 

To  derive  an  expression  for  G^  /G„   without  the  use  of 

En      Ep 

scaling    laws,    use    is   made    of   equations    previously    introduced, 
specifically 

GC   E    <%    +   GEn»    DC  (11) 

and 

F  M 

GM   e    (G         +    G_     )    D^    +  (GM      +    GM    )2    D™  (13) 

M  Ep  En '       M  y    Mp  Mn '         M 

where    D    ,    D    ,    and   D      are    defined    in    equations    (6),     (8),    and    (9) 
respectively,    and 

G2    =   A(q2)     +    B(q2)    tan2   § 
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2  2  2  2 

A(q    )     =    GCD    +    -T?DGMD 


B(^)     =   f^D^D 


(5a) 


By   the  above   equations, 
.2 


2       2 

G~    =    (G         +  G       )       D„ 
D         v    Ep  En;         C 


1    2 
—  T)     T(G         +    G       )DE    +     (G         +   G       )2DM         f(0)        (36) 
3     '     [K    Ep  En;    M         V    Mp  Mn;       Mj         v     '       K^    ' 


where 


f(0)    =1+2    tan2  §■    . 


For    the  proton,    use    equation    (26) 

_2 


g2  =  t£-  + 

P         X^p 


V. 


r-2~   f(Q)    G 


Mp 


and  recall    from    equation    (18) 


/         \ 


R    = 


/  \     OVP  f  »     OYT> 


2      2 
By    taking    the   ratio   G  ,/G      it    is   possible    to   rearrange   terms    to 
•^  dp 


find 

4 


Mott 

J? 


(G   +G   )2D2+?-?7  |(G   +G   )DE+2(G   +G   )  DM 
v  En   En ;   C  1  7  LV  En   En  '  M   v  Mn   Mn '    ! 


d0_ 

d/u 


Mott 


D 


Ep       Eny       C   3  '    Lv    Ep      En;    M       v    Mp      Mn '    M      f(Q) 

fk.  +  Jk.  f(G)  G2 

1+7?p     1+7?p  MP 


Again   rearranging    terms,    end   up   with, 
Mott 


R 


dA 


da 

dA 


Mott 


£IS1  +  Jk_  f  2me 


D2(^+l|       P(0) 


D 


2«      / JB*   /^EnY       o^2/GEnV.     E   Mf   GEr 


(37) 


/g     +g 

Mp      Mn 
V         Ep 
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G..    +G. 


E  /J    J         Mp      Mn  1       .  ,_M.2  /     Mp_Mn 


where 


-1 


F(9)    =    f(0) 

This    equation    is    free   of  approximations   and    does    not    contain 
scaling . 

To  put    equation    (37)    in    to  a    form   free   of   magnetic    contri- 

2 
butions,    a    resort    to   graphical   methods    at    constant    q      was    made. 

The   equation    can    be  written    in    the   form 


Y    =   mX   +   b 


where 


Y  =   R 


d_£ 


Mott 
P 


d£_\MOtt 


mi  +  _2 

l-Kfl  1-tT) 

p 


2 
G 

4> 


m   =    D 


En 


Ep 


+    1 


X   =    F(0)    =    (1    +   2    tan2   |) 


-1 


b   =  |-      (/   fen    )    +    2y   S*    +   ADE    pM  /°En  \  fap^ln 
3      DlM      VW  "      %  "      MiGEp/lGEp        J 

*   De2+   4DE    DM  (?^G"A   +    4DM2(/GMp+GMn 

M      M   M  V   %   /      "  \   GEp 

2 
For  constant  q  ,  X  and  Y  are  the  only  variable  terms,  Y  =  f(R,0) 

2    2 
and  X  =  f(6).   The  assumption  was  made  that  the  error  in  G.^/G^ 

v  '  ^  Mp   Ep 

2         -2 
for  0.1  £  q  ^  0.8  F   has  a  negligible  effect  when  compared  to 

2    2 
the  error  in  R .   It  should  be  noted  here  that  G   /G„   has  to  be 

Mp   Ep 

obtained  either  from  the  scaling  laws  of  the  proton,  or  f  r  on  the 

de  Vries  b'  fit;  a  later  discussion  shows  that  it  makes  little 

difference  which. 
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Using  a  least  squares  fit  to  the  data  for  the  above  equations, 
the  graphical  methods  were  found  to  be  relatively  inaccurate. 

The  remedy  to  the  above  difficulty  was  to  turn  to  an  analytical 
method  where  the  points  were  averaged  using  a  weighted  average 
method.   Equation  (37)  can  be  rearranged  to  give 


R 


/d£_f0tt 
ld"/p 

'do\Mott 


F(Q)    +  J!p_  Se 

1%     ^P  4p 


-^   )[D"    F(9)    +  -  77    D1" 


En 
'Ep 


2  h.  E         8  EM 

2DcP(6)+1r,DDM    ♦   3nDDMDMS, 


(38) 


where 


+    fj?    F(e)+   |  ,,       DE2+    8  DE      M  8  DM2    2 

C       V     ;       3     'D      M         3      D      M      M  3     'D      M 


G       +G 

s  =    "r    "n 


Ep 


The  assumption   was    made   that    every   quantity   can    be   taken    to   be 

exact    (although    the    D's    are    based   on   theoretical   models),    except 

2         2  2         2 

R,    G,,   /G^,    S,    and   consequently    the   unknown,    G^   /G^ 
Mp      Ep  ^  En       Ep 

Up    to    this   point    equation    (38)    is    exact   and   contains    no 

2         2 
scaling,    however,    G„   /G„    ,    and   S   are  as    yet    undetermined.      There 

Mp      Ep 

2         2 
are   two  methods   available    to   get    G      /G        since   absolute    data   are 

Mp   Ep 

not  available.   The  methods  are  from  the  de  Vr ies  b'  fit  and  the 
scaling  laws . 

If  the  de  Vries  b'  fit  is  used, 

Y 


%H%Uviies 


GMp+GMn 
— - % )  de   Vries    . 

GEP 
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If   scaling    is    used, 


=  p. 


S  =  M-      +  M- 
P  " 


Now   equation    (3$)    is   put    into   a   more   workable   quadratic    for 

aX2    +    bX    +    c    =   0 


m 


where   X   =   G^  /G„    .       By    inspection 
En      Ep 


a    =    D2C  F(9)    +   f  nD  < 


b   =    2DcF<e>+!%DM+f"DDM    °M    S 


2  .        ^E2       8  JS    31    „         8  M2       2 


c    -    Dc    F(0)     +   -  f|D   DM    +    -  r?D    DM    DM    S    +    3  7]D    DM       S 


Mott 


rdQ 


Mott 

! 

D 


r? 


l+n       v    ;       i+ti 

p  p 


GEp, 


The   results    of   these    two  methods    are    summarized    in   Table   VII. 

By   applying    Schumacher's    scaling    to   equation    (38),    a    comparison 

of   G„   /G„      can   be  made    to   that    of    the   two   methods    discussed  above. 
En      Ep 

/  ^  o  \                                  •          -■■         ~2     .2         „           Mp      Mn 
In   equation    (3°),    the    terms    involving   G      /G       ,    S   =  — 'p. , 


Mp   Ep 

2 
and  S   are  the  only  ones  affected.   Let  X  =  G„  /G 

En   Ep 


application  of  equation  (4),  it  is  seen  that 

2 


5^V 


=  (p.  +  H  X) 
v  p   n 


and 


S  =  (Up+  M-n)(l  +  X)  . 


Ep 
Then  by 


(39) 
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By   rearranging    terms    after    substituting    equation    (39)    into 
equation    (3&),    "the   resulting    equation    is 


wher  e 


2  2  E 

a*    =    D      F(9)     +  -  7?       D 
C      K     '  3      D      M 


a 'X     +    b!X   +   c'    =   0 
2 


(40) 


-R 


dsi 


Mott 


r? 


d£_ 


Mott 


D 


2       8  EM    .  .8  M    .  .2       //A    . 


1+ri         n       3       D  M   M    v    p      n'       3      D   M    v    p      n 
P 


b'    =    2D2C   F(0)+   |^7D   D*    +2 


8  EM 

—  r?      D     D      (u.   +  u    ) 

3      D     M      M    y    p         n7 


+    2 


and 


8  ^M    /  x 2 

3      D      M    x    p      n' 


d£ 
-R iE 


Mott 


:!e_ 


p,    m- 


vMott       Htjp        p    n 
dsi  J 


C*    =c    =    nF(9)+—  T7      D      +  —  r)      D      D'      (\i      +  u    ) 
C       v     ;  3      D      M         3      D     M      M    v    p  n; 


8  M  2 

3      D      M         p  n' 


(40b) 


da 


Mott 


„  1  dA  / 


dCJ 
dA 


Mott 


1*1?        m?p      p 


(40c) 
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III.       REEVALUATION    OF    EXPERIMENTAL    DATA 

A.  DATA 

The  data  analyzed  here  were  collected  by  Burailler  and  others  Ll9J 

at  the  Naval  Postgraduate  School  linear  accelerator  in  the  range 

2         -2 
0.1  £  q  ^  0.4  F   ,  and  at  the  Mark  III  accelerator  at  Stanford 

2         -2 
University  in  the  range  0.2  s,  q  <,    0.8  F   .   The  experimental  set 

up  was  basically  the  same  in  both  laboratories.   Whenever  possible, 
overlapping  sets  of  data  were  taken  to  check  for  any  systematic 
errors  in  the  equipment.   No  systematic  deviations  between  the  two 
sets  of  measurements  were  found  and  the  data  analysis  was  per- 
formed independent  of  the  source.   Table  I  shows  the  data  used  in 
these  analyses,  with  the  error  listed. 

B.  ROSENBLUTH  PLOTS 

Equation  (17)  formed  the  basis  for  the  Rosenbluth  plots.   This 

equation  was  broken  up  as  follows  for  the  plots: 

theor     (  ,  „  \Mott 

'£•)„ 


m: 


2    0 
X    =    tan      - 


where 


cer- ■  «*-.  Mr 

and   R    is    given   by    equation    (18).      Wherever    appropriate,    weighted 
averages   were   used,    and   the   corresponding    error    calculated.       A 
weighted   least    squares    fit   was    used   to    determine    the  slope   and 
intercept,    with   the   error    in    each.      This    method    is    described    in 
Appendix   C. 
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2 
Since   these  plots    were   made   for    a    constant    q    ,    it    is    readily 

apparent    that    from   the   slope  m, 

°MD   (q    }    "    4  77  D 
and   from   the   slope  ra  and    intercept    b, 

GED    (q    )    =    b    -    %m. 

Table  III  gives  the  results  of  these  plots  for  the  deuteron, 
and  a  sample  plot  is  shown  in  Figure  1. 

C.  MAGNETIC  CORRECTIONS 

The  magnetic  correction  gives  the  number  by  which  the  experi- 

2      2 
mental  G   or  G  ,  is  to  be  multiplied  to  remove  the  magnetic  contri- 
p      d  ^ 

2       2 
butions  and  yield  G„   or  G_  ..   These  corrections  were  computed 

*  Ep  Ed  K 

using    equations    (22),    (25),     (27),    (32),     (34),    and    (35).      The    Dc 

terra   contained    in    equation    (32)    is    the   charge   structure   factor 

and    is    given   by    equation    (6),    with   the   Feshbach-Lomora  wave 

functions    LllJ    being   used   to    describe    the   deuteron.       Values    of    D 

^  c 

are  listed  in  Table  III.   The  magnetic  corrections  using  both 
scaling  laws,  a  comparison  of  the  resulting  corrections,  and  the 

effect  of  magnetic  corrections  using  Schumacher  scaling  on  G^  /G„ 

v  v  v  En      Ep 

as    calculated  by    Stewart   are    shown    in   Table   VI.      The   error    was 
calculated   by    taking  K   /K      as    the  standard   value. 

D.  CALCULATION  OF  G^  /G^      WITHOUT  SCALING 

En   Ep 

Equation  (38)  is  the  basis  for  these  calculations.   The  D  , 

E       '  M 
D  ,  and  D   are  respectively,  the  charge  structure  factor,  magnetic 

contribution  to  the  scattering  process  arising  from  the  convection 

of  charge  in  the  deuteron,  and  the  contribution  of  the  intrinsic 
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magnetic   moments    of   the  proton   and   the   neutron    to   the   scattering 
process.      These   factors    are   given   by    equations    (5),    (8),    and    (9). 
As    with   the  magnetic    corrections,    the    Feshbach-Lomom  wave   functions 

were   used.      The    values    of    these   wave    functions    at    the  appropriate 

2 
q      are  also   listed   in   Table    III. 

2         2 
G      /G„      is    contained  within    equation    (38)    and    is    an   undetermined 
Mp      Ep 

quantity.      Since  absolute   data    are   not   available,    this    ratio  may  be 

found   from   the    de  Vries   b'    fit,    or    from   the  proton    scaling    law. 

Both   methods    of    calculation   were   used,    and    it   was    found   that    they 

differ    by   a    small   percentage. 

Results    of   these    calculations    are   listed    in   Table   VII,    and   a 

comparison    of   G^   /G^      computer    using    data    from   de   Vries    b'    fit, 
En      Ep 

and  using   proton   scaling    is    shown    in    Figure    5. 


E.       CALCULATION   OF  G^    /G_      USING    SCHUMACHER    SCALING 

En   Ep 

These  calculations  were  made  using  equation  (4-0),  with  the 

coefficients  given  by  equations  (40a-c).   As  in  the  previous  cal- 

E       M 
culations  the  D  ,  D  ,  and  D   were  calculated  using  the  Feshbach- 
c   M       M 

Lomom  wave  functions  as  listed  in  Table  III. 
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IV.   PRESENTATION  OF  RESULTS,  DISCUSSION,  AND  CONCLUSIONS 

A.   PRESENTATION  OF  RESULTS,  AND  DISCUSSION 

The  experimental  data  are  presented  in  Table  I.   The  errors 
quoted  are  attributed  to  counting  and  target  thickness  errors, 
and  are  the  only  errors  that  enter  into  a  ratio  experiment. 

Table  II  shows  the  Mott  cross  sections  for  point  proton  and 

2   2 

deuteron  in  the  laboratory  system;  G  (q  )  from  the  b'  fit  of 

de  Vries  is  the  absolute  form  factor  for  the  proton;  the  theo- 
retical cross  sections  as  described  in  Section  IIA;  and  the  nor- 
malization factor  N  required  to  make  the  experimental  deuteron 

cross  section  absolute  in  the  first  order  sense. 

E    M 
Shown  in  Table  III  are  the  quantities  D  ,  D...  Dw  as  described 

^  c        M        M 

in    Sections    IIA,    computed   from   the    Feshbach-Lomom  wave    functions 

without   relativistic    corrections. 

Table    IV   lists    the  magnetic   and   electric   form   factors    of    the 

deuteron   as    calculated   from    the  Rosenbluth   plots.       Figures    2    and 

3    illustrate   these   results.       It    is    seen    from   Figure    2    that    the 

electric    form   factors   approach    the    static    limit    as   predicted   by 

the   theory.      From   Figure   3    it    is    seen    that    the  magnetic    form 

factors    are   to   large   by   a    factor    of   three.      This    is    due    in   part 

to   the    fact    that    G        is   a    function    of    the   slope    of    the   linear    plot 

MD  v  r 

and    is    very    senstive   to    experimental    errors.      A  weighted   least 

squares   plot   was    used   in   an    attempt    to   minimize    the    error.      Also 

G 

shown  in  Table  IV  is  the  ratio  —■  ■■  —  .   If  scaling  is  to  hold, 

u,  '  G 
D  ED 

this  ratio  must  be  unity,  within  experimental  error.   It  is  sepn 
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MD 
however,  that  within  experimental  error  — — ^  3,  indicating 

^D  GED 

that  scaling  does  not  hold.   Since  the  values  of  G    are  about 

MD 

three  times  larger  than  predicted,  it  must  be  held  that  this  test 
is  inconclusive. 

To  investigate  the  effect  of  magnetic  contributions  to  the 

electric  form  factor,  and  get  information  as  to  the  importance 

2 
of  scaling  within  the  range  of  q   of  interest,  the  Rosenbluth 

plots  were  considered.   From  these  plots  it  was  seen  that 

4d  =  f<b-,r-> 

where  b    is    the   Y   intercept,    and   m    is    the   slope   of  the   linear    plot. 

Specifically, 

GED  =  b    -  h  * 

where  m   contains    the  magnetic    contribution.       Define 

*2 

G™    =    b 

ED 

■* 
where  G    is  the  electric  form  factor  without  magnetic  contri- 
ED  y 

butions.   Table  V  shows  the  results  of  these  calculations.   The 
result  is  that  a  very  small  error  is  introduced  if  the  magnetic 

contributions  are  disregarded,  thus  concluding  that  scaling  has 

2 
little  effect  in  the  range  of  q   considered  here. 

Stewart  LlO]  was  interested  in  the  ratio  G  „/G„  ,  from  which 

ED   Ep 

he  finds  G^   /G^_ .       It  was  shown  in  Section  IIB,  that  the  ratio 
En   Ep 

G-^/G^   contains  the  magnetic  correction  ratio,  which  in  turn 
ED   Ep 

contains  the  scaling  laws.   Using  both  scaling  laws,  a  comparison 
was  made  of  the  two  magnetic  correction  ratios.   Table  IV  contains 

these  results,  showing  that  the  difference  in  the  corrections  is 

2 
small,  but  tend  to  increase  with  q  . 

Table  VII  shows  four  sets  of  values  of  G„  /G„   as  obtained: 

En   Ep 

(a)  Without  using  scaling,  but  using  the  assumption  of  the  de  Vries 
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b'    fit:     (b)    By    using  proton    scaling    in    the   form  G        =  M-   G      : 
K     '  Mp  p    Ep 

(c)    By   Using    Schumacher's    scaling    law   for    the  proton   and;     (d)    As 
calculated  by    Stewart   LlQj  .       In    each    case   the   Fes hbach -Lorn ora  wave 
functions    were   used,    and  relativistic    corrections    described  by 

Stewart    were  applied.      Each    of   the    columns    shows    the   slopes 

2 
d(G      /G     \/dq    ,    calculated  by    using   a    weighted  least    squares    fit. 

In    each    case   the   slopes    are    in    good  agreement. 

B.       CONCLUSIONS 

Based  upon  the  experimental  data  available  in  the  range  of 
momentum  transfers  of  interest  here,  the  following  conclusions 
can  be  made: 

(a)  Both  forms  of  scaling  as  applied  to  the  deuteron  form 

factor  have  negigible  effect  on  the  ratio  G^  /G^  . 

En   Ep 

(b)  The  application  of  either  scaling  law  to  the  proton  form 

factor  results  in  nearly  identical  corrections  to  the  ratio 

G„  /G„  ,  thus  allowing  no  distinction  between  them. 
En   Ep 

(c)  The  agreement  of  the  experimental  data  with  the  neutron- 

2 
electron  interaction  slope  at  q  ~  0  is  not  influenced  by  the 

choice  of  the  scaling  law. 

Although  both  forms  of  scaling  give  satisfactory  results  in 

the  momentum  transfer  range  investigated  here,  the  Schumacher 

scaling  appears  to  be  on  a  better  theoretical  basis,  and  should 

be  tested  at  higher  momentum  transfers.   A  conclusive  test  for 

the  scaling  laws  would  require  an  improvement  in  experimental 

accuracy  by  at  least  an  order  of  magnitude,  a  technically  nearly 

impossible  requirement. 
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The    interpretation    of   measurements    at    higher    momentum   transfers, 
where    the   difference  between    the   scaling    laws    would  be  more  pro- 
nounced   is    hampered  by    increased   difficulties    in    the   theoretical 
prediction    of    deuteron  wave    functions. 

In   agreement   with   Stewart,    an    important    qualitative   conclusion 
from   these   analyses    is    that    the   neutron    has    a    non-zero    charge    form 
factor,    not   affected  by    the   choice    of   scaling    laws,    implying   a 
charge    distribution   within   the   neutron.      The   neutron   will    then   have 
a    charge   structure  with   the    outermost  part    of   the   charge   distri- 
bution  being    negative    in   sign. 
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APPENDIX  A 
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APPENDIX  C.  STATISTICAL  TREATMENT  OF  DATA 

A.   WEIGHTED  METHOD  OF  LEAST  SQUARES 

The  principle  of  least  squares  can  be  stated  as  follows:   The 
most  probable  value  of  a  quantity  is  obtained  from  a  set  of  mea- 
surements by  choosing  the  value  which  will  minimize  the  sum  of 
the  squares  of  the  deviations  of  these  measurements.   Specifically, 
if  two  variables  X  and  Y  are  known  to  be  related  by  a  linear 
equation  of  the  form  Y  =  mX  +  b,  where  m  is  the  slope  of  the  line 
and  b  its  Y  intercept,  and  if  a  series  of  N  observations  (x.,y.) 
are  made  in  which  random  errors  occur  only  in  the  y.  measurements, 
and  the  observations  are  not  all  samples  of  the  same  parent  dis- 
tribution, then  the  deviations  are  weighted  inversly  as  the 
variances  of  their  parent  distributions,  and  Young  L17J  shows 
that  the  quantity  to  minimize  is 

2      4 

£  w.d.  =  S  — —  =  T,    w.(rnx.+b-y  .)  (4l) 

11       2        l v   l     iy  v   ' 

S. 

l 

where 

w.  a   — - 

1   s2 

yi 


and  S         is    the   variance    in    y.    Ll8J.       Upon   minimizing    equation    (4l) 
i 

with      respect    to   b  and  m,    the   normal    equations    are 


mSw.x.    +bSw.    =  S   w . y . 
li  l  li 


2 
mSw.x.    +   bS    w.x.    =£   w.x.y.     . 
li  li  ill 


Solving   for    m   and   b,    it    is    found    that 
m    =   M/D 


(42) 


(43) 
b    =    B/D 
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where 


M   = 


E   w.y.  T,   w . 

11  1 


S    w.x.y.        T,    w.x. 
111  11 


B 


S    w.x . 
i    i 

E   w.x. 
l    l 


£   w.y. 
l    l 


T,   w.x.y. 
ill 


and 


D 


S    w.x . 
i    i 

^  2 

ii    w.x. 

1    1 


£   w 


£   w.x. 
i    l 


The  variances  in  m  and  b  can  be  calculated  by  use  of  the 
equation 


sq  =  e(& 


(44) 


2  2 

where  S^  is  the  variance  of  the  mean  of  Q,  S    is  the  variance 

Q  y± 

of   the  mean    of    y.,    and   so   forth.      Application    of   equation    (44)    to 
the  present    situation    gives 


and 


where 


£    w.x. 
^2         ~2  11 

S,     =    S        — 

b  y  D 


2  2      S    Wi 

S      =    S  _ 

m  y  D 


c^2  r^2 

S      =   w.    S 

y        i    y. 


(45) 


B.       WEIGHTED  AVERAGES 

If   a    quantity   x    is   measured  by   N   methods,    obtaining    values 

X.  .  .  .x    ,    which   are   known    to   have   standard   deviations    S,...S.T, 
IN  IN 

then   a    simple   arithmetic   average  would   not    give   the   best    value 
because    it   would  make   the   x.    equal    in   weight,    while    their    errors 
are  assumed   to  be    different.      A   weighted  average    is    then   used  as 


44 


a  means  of  determining  the  best  value.  The  question  is  how  to 
assign  the  correct  value  of  the  weight  factor  w.  to  x.  to  give 
the  minimum  standard  deviation  in  the  average  value  x.  Making 
a  generalization  of  Beers  Ll8]  it  is  seen  that 

£  w.x . 


2  w . 

1 


and 


wher  e 


sf 

= 

1 

^       1 

X 

X 

w. 

1 

0! 

1 

s2 

1 

C.       PROPAGATION   OF   ERROR 

Young   Ll7l    shows    that    if   a    quantity    Q   is    to   be    calculated 
from   several    observed   quantities    x.    and 

Q  =  f(x±) 

then 

*0-*($)\  (47) 

2  2 

where  S^  is  the  variance  of  the  mean  of  Q,  and  S    is  the 
Q  x.  . 

variance  of  the  mean  of  x..   Equation  (47)  can  be  shown  to  be 

true  even  if  different  numbers  of  observation  are  made  on  the  x. 

l 

Application  of  equation  (47)  gives  the  variance  of  the  mean 
for  the  appropriate  functional  form  of  Q. 
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